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We discuss the excess in the hadroproduction of B mesons at the Tevatron. We show that an 
accurate use of up-to-date information on the B fragmentation function reduces the observed excess 
to an acceptable level. Possible implications for experimental results reporting bottom quark cross 
sections, also showing an excess with respect to next-to-leading order theoretical predictions, are 
discussed. 



Since a few years, bottom production has been one 
of the very few instances in which experimental results 
and Quantum Chromodynamics (QCD) predictions have 
sometimes displayed not too good an agreement. Bot- 
tom quark hadroproduction cross sections have been 




at the CERN 
and DO il 



measured by the UAl Collaboration 
SppS collider and by both the CDF [| 
experiments at the Fermilab Tevatron in pp collisions, 
and found to be about a factor of two or more larger 
than next-to-leading order (NLO) QCD predictions 
CDF has recently published data for i?+ meson produc- 
tion They claim an excess over QCD predictions by 
about a factor of three. The HI and ZEUS experiments 
at the electron-proton collider HERA have both mea- 
sured D* production cross sections pd] , p^ . Both mea- 
surements are compatible with QCD calculations [ p^ , 
although the ZEUS data is on the high side of the the- 
oretical uncertainty band. Bottom production has been 
measured at HERA and from photon-photon collisions at 
LEP, and found to be larger than predictions, by about a 
factor of three or more ||l4[--[17| . The cross section for the 
production of the heaviest quark, the top, has also been 
measured at the Tevatron, and found instead perfectly 
compatible with theoretical expectations. 

By pushing the parameters of the theoretical calcu- 
lation to somewhat extreme values, it is not impossible 
to accommodate the bottom spectrum observed at the 
Tevatron. Alternatively, one can take the discrepancy 
more seriously, and invoke some "new physics" contribu- 
tion [|l8| in order to explain it. It has also been known 
for a long time that the fixed order, NLO QCD calcula- 
tion may be insufficient to explain the data, because of 
the presence of some enhanced contributions, that can 
be included via resummation of large classes of Feynman 
diagrams, and that contribute positively to the cross sec- 
tions. These contributions are threshold effects, small- 
X effects, and high transverse momentum logarithms, 
which may be important since most of the cross section 



is measured at large transverse momentum. A full cal- 
culation of next-to-next-to-leading QCD contributions, 
years ahead in the future, might finally also contribute 
to explain the apparent discrepancy. 

In this Letter we shall not try to improve on the per- 
turbative aspects of heavy quark production. We shall 
instead focus our attention on a specific non-perturbative 
issue, namely on the implementation of hadronization ef- 
fects. In fact, we shall argue that a good part of the dis- 
crepancy between theory and data arises when one tries 
to supplement the perturbative prediction for h quark 
production with a non-perturbative model for the forma- 
tion oi a B meson from the h quark, or, alternatively, to 
correct the data in an attempt to give a h quark spec- 
trum rather than a B meson one. The non-perturbative 
hadron formation effect is usually introduced by writing 
the hadron-level cross section for B mesons as 



— — — [ dpT^dz— — D(z) 6{pt — zprj.) 
dp-v J dpT 
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the function D{z) being a phenomcnological parametriza- 
tion of hadronization effects. Traditionally, the Peterson 
et al. |l9) D{z; e) form of the fragmentation function is 
used, implemented in conjunction with a quark cross sec- 
tion given by a shower Monte Carlo program. The e pa- 
rameter is obtained from fits to e+e" data The 
effect of fragmentation is to reduce the momentum of 
the B meson with respect to that of the b quark. It is 
roughly a 10% effect, being of the order of A/m, where A 
is a hadronic scale, of the order of a few hundred MeV, 
and m is the bottom quark mass. It has however an im- 
portant impact on the value of the cross section, because 
of the steeply falling transverse momentum spectrum of 
the b quark. Since transverse momentum cuts are always 
applied, the measurable cross section is strongly reduced 
by this effect. It should be clear from this discussion that, 
in order to assess the presence of a discrepancy in the B 
production data, the effect of fragmentation should be 
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assessed clearly and unambiguously. 

In Ref. |lO| the CDF Collaboration compares its data 
to a theoretical prediction obtained by convoluting the 
NLO cross section for bottom quarks with a Peterson 
fragmentation function. They use e = 0.006 ± 0.002, 
which is the traditional value proposed in Ref. |Q . They 
claim that their data is a factor of 2.9 higher than the 
QCD calculation. 

The purpose of this Letter is precisely to implement 
correctly the effect of heavy quark fragmentation in the 
QCD calculation. Several ingredients are necessary in 
order to do this: 

• A calculation with resummation of large transverse 
momentum logarithms at the next-to-leading level 
(NLL) should be used for heavy quark production 

, in order to correctly account for scaling viola- 
tion in the fragmentation function. 

• A formalism for merging the NLL resummed results 
with the NLO fixed order calculation (FO) should 
be used, in order to account properly for mass ef- 
fects 



in the following. 

• A NLL formalism should be used to extract the 
non-perturbative fragmentation effects from e~^e~ 
data l2Mgtl. 



We begin by pointing out that, as shown in Rcfs. P7| , p8| , 
the value e = 0.006 is appropriate only when a leading-log 
(LL) calculation of the spectrum is used, as is the case in 
shower Monte Carlo programs. When NLL calculations 
are used, smaller values of e are needed to fit the data. 
It must further be pointed out that, as noted in pO| , ^ , 
it is not the detailed knowledge of the whole spectrum 
of D{z) in z g [0, 1] to be relevant for the calculation of 
hadronic cross sections. For the steeply falling differen- 
tial distributions da/dp-n that have usually a power law 
behaviour, the knowledge of some specific moment of the 
fragmentation function 



(2) 



is sufficient to obtain the hadronic cross section. In fact, 
assuming that da/dp^: — ApT^"^ in the neighborhood of 
some value, one immediately finds 



da f , , /\^r/ ^T-, 
- — = / dzdpT D[z)— d{pT - zp-r) = —Dr, 

dp-T J Pt Pt 



(3) 



Thus, the hadronic cross section is given by the product 
of the partonic cross section times the n*^ moment of the 
fragmentation function, where n is the power behaviour 
of the cross section in the neighborhood of the value of 
being considered. In Ref. |3l] it is also shown that this 
is an excellent approximation to the exact integral in the 
cases of interest. The value of n for the p-r spectrum in 



the region of interest ranges from 3 to 5. It is therefore 
clear that, when fitting e+e~ data, getting a good deter- 
mination of the moments of the non-perturbative frag- 
mentation function between 3 and 5 is more important 
than attempting to describe the whole z spectrum. 
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Moments of the measured B meson fragmentation 
function, compared with the perturbative NLL calculation 
supplemented with different D{z) non-perturbative fragmen- 
tation forms. The solid line is obtained using a one-parameter 
form fitted to the second moment. 

Fig. |l| shows the moments calculated from the xe (the 
B meson energy fraction with respect to the beam en- 
ergy) distribution data for weakly decaying B mesons 
in e+e^ collisions published by the ALEPH Collabora- 
tion . The experimental error bars shown in the plot 
have been evaluated by taking into account the full bin- 
to-bin correlation matrix p3[ |. Four curves are superim- 
posed to the data. All of them have been obtained with 
an underlying NLL perturbative description j2^j2^ . The 
bottom quark mass m has been taken equal to 4.75 GeV 
and the QCD scale has been fixed to A^^) = 0.226 GeV. 
Sudakov resummation has not been included, since its 
effect is negligible in the low-moment region |2^. These 
arc the default values of the parameters that we shall use 
in this work for the computation of the hadronic cross 
section. 

The dot-dashed line represents the purely perturbative 
part. The dashed line represents the convolution of the 
perturbative part described above with a Peterson form 
with e — 0.006. It is evident that this produces a poor 
description of even the lowest moments. The dotted line 
is obtained using e = 0.002, a value known to produce 
good fits of the xe distribution when used together with 
a NLL perturbative calculation |^,^. The description 
of the moments improves, but the line still cannot fall 
within the error bars. There is thus a problem in obtain- 
ing a good fit of the low moments of the fragmentation 
function using the Peterson parametrization. The prob- 
lem can be traced back to the need to fit points with very 
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large xe (where most of the e+e" data is) since there the 
perturbative calculation becomes less reliable. Normally, 
the very large x e region is excluded from the fit because 
of this reason. The computed cross section is thus al- 
lowed to become negative in this region, a fact that leads 
to an underestimate of the low moments. 

It should be clear from the aforementioned arguments 
that, in order to make accurate predictions for hadronic 
cross sections, the non-perturbative part of the fragmen- 
tation should be fitted in such a way that the low mo- 
ments are well reproduced. This is shown in the solid 
line in the figure. A one-parameter form of the non- 
perturbative fragmentation function has been used^ and 
its free parameter has been fixed by fitting the N = 2 
point in moments space, i.e. the average energy fraction 
(xe)- In this case the functional form is good enough to 
describe well the experimental data up to ~ 10. It is 
therefore a good candidate to be employed in the calcu- 
lation of the hadronic cross section according to Eq. (|^). 
We shall refer to this fit in the following as the "N — 2 
fit". 

We notice that the effect of non-perturbative fragmen- 
tation (i.e., the ratio of the dashed, dotted and solid 
curves with the dot-dashed curve) is considerably re- 
duced when the moments are fitted. Thus, perturbation 
theory alone gives a much better description of the low 
moments of the fragmentation function, rather than of its 
shape in x-space, requiring less non-perturbative input. 
This is a consequence of the fact that, at large x, many 
enhanced non-perturbative contributions and hadroniza- 
tion effects due to the limited phase space come into play. 
Indeed, the form of the leading power correction in mo- 



ments space is well known, and reads |35 
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FIG. 2. Prediction for the B cross section, obtained using 
the calculation of Ref. |Q supplemented with the A'' = 2 fit 
of the non-perturbative fragmentation function, compared to 
the CDF data of Ref. For comparison, the result ob- 

tained using a Peterson form with e — 0.006 is also shown. 
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FIG. 3. Data over Theory ratio for B production. Data 
points and theoretical curves are as in Fig. ^ 



It can easily be checked that the form we employed in the 
A^ = 2 fit is consistent with this leading power correc- 
tion, provided one replaces a with 2m/ A. One can also 
clearly see the non-perturbative correction to be mini- 
mal for N = 2. It is therefore always desirable to study 
moments, rather than the x shape of the fragmentation 
function, also in order to perform QCD studies. 



In Figs. H and ^ we show the final prediction for B 
hadroproduction at the Tevatron, obtained by the proce- 
dure outlined above. It is clearly shown how the "Peter- 
son with e = 0.006" choice underestimates the B cross 
section at large values of Pt- It is also clear that the 
claimed discrepancy of a factor of 2.9 is now reduced to 
a factor of 1.7 with respect to the central value predic- 
tion. The band obtained by varying the scales gives an 
idea of the theoretical error involved, and, as one can see, 



^Exactly which functional form is used is actually not rele- 
vant. However, for the sake of completeness, we mention that 
we have used a normalized Kartvelishvili et al. form | p^ , 
D{z- a) = ia+ 1)(q + 2)z"(l - z). 
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the data are not far above it. 
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FIG. 4. The effect of the different ingredients in the cal- 
culation presented in this work, relative to a fixed order cal- 
culation with Peterson fragmentation and e = 0.006. 

At this point, we wish to quantify to what extent the 
various ingredients of the present calculation affect the 
computed cross section, so that it is in fact larger than 
the one given in Ref. . This is shown in Fig. ^. From 
the figure we see that the FO (Fixed Order) calculation 
with e = 0.006 is the lowest curve. Using the more ap- 
propriate value e = 0.002 brings about a 20% increase of 
the cross section at Pt = 20. Using the FONLL calcula- 
tion of Ref. brings about another 20% increase, and 
so does also the use of the = 2 fit. The total effect 
is an increase by a factor of 1.2^ ~ 1.7, which turns the 
factor of 2.9 reported in Ref. into the 1.7 observed 
here. 

Our FO, e — 0.006 cross section is also higher than 
the one presented in Ref. llo[| in the low region. This 
difference could be due to the different possible treat- 
ments of fragmentation at small transverse momentum. 
We have applied the fragmentation to the momentum, 
rather than the energy or the -f component, of the frag- 
menting particle. We believe that these other choices, 
although acceptable in the large-p^ region, are not ap- 
propriate in the non-relativistic limit. 

The SLD experiment has also published accurate data 
on the h fragmentation function . Using their data in- 
stead of the ALEPH ones brings about a slight decrease 
of the cross section, below 4% in the region of interest. 
On the other hand, using more recent parton distribution 
function sets [3^ ], we find an increase of the predicted 
cross section between 4 and 8% in the region of interest. 

In summary, we find that an appropriate treatment of 
the fragmentation properties of the h quark considerably 
reduces the discrepancy of the CDF transverse momen- 
tum spectrum for the B mesons and the corresponding 
QCD calculation. The experimental points are compat- 
ible with predictions obtained using the present value 



of the QCD scale parameter and of the structure func- 
tions, and a 6 pole mass of 4.75 GeV, and lie near the 
upper region of the theoretical band obtained by varying 
the factorization and renormalization scales. Including 
experimental and theoretical uncertainties, the updated 
Data/Theory ratio can be written as 1.7 ± 0.5 (expt) ± 
0.5 (theory). The calculation we have adopted includes 
in a consistent way fixed order QCD results and the NLL 
resummation of transverse momentum logarithms. Fur- 
thermore, the "moments" method introduced here avoids 
the difficult large-x region in the fragmentation function, 
that would require more complex treatment and intro- 
duce further uncertainties. 

While we have here convoluted a perturbative pre- 
diction to get a hadron level result, the opposite path 
is followed by experiments when they deconvolute their 
hadron level cross sections in order to publish quark level 
data [|l]-||. In the light of what argued in this Letter, and 
of the apparent excess also shown by those data, it will 
be advisable to investigate whether a similar bias might 
have affected those results. 

In the meantime, we emphasize the importance of the 
direct measurements of the moments of the fragmenta- 
tion function for heavy quarks. Measuring directly mo- 
ments instead of the x distribution could be useful for the 
purpose of QCD studies, and also for the computation of 
production cross sections in hadronic collisions. 
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